The Simulation of the Inelastic Impact 

Hiroto Kuninaka* and Hisao Hayakawa 1 

* Graduate School of Human and Environmental Studies, Kyoto University, Sakyo-ku, Kyoto, Japan, 606-8501 
' Department of Physics, Yoshida-south campus, Kyoto University, Sakyo-ku, Kyoto, Japan, 606-8501 

Abstract. The coefficient of normal restitution (COR) in an oblique impact is theoretically studied. Using a two-dimensional 
lattice models for an elastic disk and an elastic wall, we investigate the dependency of COR on an incident angle and 
demonstrate that COR can exceed one and have a peak against an incident angle in our simulation. Finally, we explain 
these phenomena based upon the phenomenological theory of elasticity. 



INTRODUCTION 



The coefficient of restitution (COR) e is introduced to 
determine the post-collisional velocity in the normal col- 
lision of two materials and defined by u • n = — eu • n, 
where u and u are respectively the velocity of the con- 
tact point of two colliding materials before and after 
the impact, and n is the normal unit vector of the tan- 
gential plane of them. In the oblique impact, it has be- 
come clear that COR depends on the incident angle as 
well as on the impact velocity. Louge and Adams 1 1] re- 
cently reported that COR increases as a linear function 
of the tangent of the incident angle in the oblique im- 
pact of a hard aluminum oxide sphere on a thick elasto- 
plastic plate. They also suggested that COR can exceed 
1 in most grazing impacts. In this proceeding, we carry 
out the two-dimensional simulation of the oblique impact 
and investigate the dependency of COR on the incident 
angle based upon the theory of elasticity. 



MODEL 

Our numerical model is a two-dimensional model which 
is composed of an elastic disk and an elastic wallj^]. 
Each of them is composed of randomly placed mass par- 
ticles connected by nonlinear springs each other. Num- 
bers of mass particles are 400 for the disk and 2000 for 
the wall. The width and the height of the wall are %R and 
2R respectively where R is the radius of the disk. Both 
sides and bottom of the wall are fixed. Spring potential is 



described as V(x) = jk a x + jkbX , where x is a stretch 
from natural length. The values of k a are 1.0 x mc 2 /R 2 
for the disk and k a = 1.0 X I0~ 2 mc 2 /R 2 for the wall, 
where c is the one-dimensional sound of velocity. And 
we adopt kf, = k u x 10~ 3 //? 2 for each of them. Poisson's 
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FIGURE 1. 

lattice. 



The elastic disk and wall consisted of random 



ratio v can be evaluated from the strains of the band of 
random lattice in vertical and horizontal directions to the 
applied force. We obtain Poisson's ratio v of the model 
as V= (7.50 ±0.11) x 10~ 2 

In our simulation, we define the incident angle y by 
the angle between the normal vector of the wall and the 
initial velocity vector of the disk(Fig.Q). We fix the 
initial colliding velocity of the disk as |v(0)| = 0.1c to 
control the normal and tangential components of the ini- 
tial colliding velocity as v f (0) = |v(0)| sin 7 and v„(0) = 
|v(0)|cosy, respectively. From the normal components 
of the contact point velocities before and after the col- 
lision, we calculate COR for each 7. We use the fourth 
order symplectic numerical method for the numerical 
scheme of integration with the time step Af = 10~ 3 /?/c. 



RESULTS AND DISCUSSION 

Figure |2 is the COR against the tangent of the incident 
angle ^\ = tan 7 in our simulation. The cross points are 



the average and the error bars are the standard deviation 
of 100 samples for each incident angle. This result shows 
that the COR increases as *Fi increases to exceed 1, and 
has a peak around = 5.0. This behavior is contrast to 
that in the experiment by Louge and Adams 1 1 ]. 

To explain this result, we consider the correction of 
COR by the local deformation of the wall. We assume 
that the normal unit vector n to the surface of the wall 
rotates toward the incoming disk by an angle a to be- 
come n". If we define e = — (u • n)/(u • n) and e a = 
— (u •n a )/(u-n a ), the relation between e and e a be- 
comes 



e= (e a +¥£tana)/(l -I'ftana), 



(1) 



where 



-tana)/(l +^itana) and = 
(¥i -tana)/(l+¥i tana) -3(1 + e a )(^ + tana)/(l - 
/Xtana). As for > we use tne phenomenological the- 
ory for the oblique impact by Walton and Braun|3j. The 
correction angle a can be estimated by the theory of 
elasticity. If we express the contact area by a parab- 
ora, tana equals to \x c — x a \(l — 2d)/R{2 — 29) with 
9 = (l/7r)arctan(l -2v)/(jli(2-2v)), where pL is the 
coefficient of friction and x c and x a are the x coordinates 
of both ends of the contact area. jj. can be calculated from 
the simulation data through the definition jj, — \ J t \/\J n \, 
where J n and J t are the normal and tangential compo- 
nents of the impulse. The relation between jj. and *J/i are 
cross points in Fig. [3] The solid curve in Fig.|2]is Eq. Q 
with e a = 0.95 which is COR in the normal impact. Our 
numerical results can be reproduced by our phenomeno- 
logical theory. 

The relation between ji and can be explained as 
follows. We assume that jags are uniformly placed on 
the surface of the wall with the density p per unit length 
and the tangential velocity of the disk is decreased by rj 
when the disk interacts with the one jag. The tangential 
and normal impulses can be calculated by calculating the 
numberof jags the disk interacts during collision time|4] 
as /, = -rarjp | v(0) | sin y%{R/c) y/ln(4c/ \ v(0) | cos 7) 
and J n = —m(e + 1 ) |v(0) cos 7. We also assume that the 
tangential impulse decreases by J t = — mQ\(0)\ siny 
which is proportional to the initial tangential velocity 
with the proportionality constant £. Thus, jx can be 
calculated by the ratio of \J, — J, \ to /„ as 



^tan7— T7ptan77r 




4c 



v(0)|cos7 



(2) 

The solid curve in Figgis Eq.(|3J with £ = 0.317 and 
Tjp = 0.0416c/R. Our numerical results can be well re- 
produced by our phenomenological theory of the coeffi- 
cient of friction. 




FIGURE 2. Numerical and theoretical results of the relation 
between S?i and COR. 
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FIGURE 3. Numerical and theoretical results of the relation 
between 1*1 and )l. 



SUMMARY 

In the present study, we have shown that COR can exceed 
1 also in our two-dimensional simulation and depends on 
ji in the oblique impact. Our results can be explained by 
our simple phenomenological theory. 
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